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After introducing essential, qualitative concepts and results, we discuss the application 
of Dyson- Schwinger equations to QCD at finite T and fi. We summarise the calculation 
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of the critical exponents of two-light-flavour QCD using the chiral and thermal suscepti- 
■ bilities; and an algebraic model that elucidates the origin of an anticorrelation between 
the fi- and T-dependence of a range of meson properties. That model also provides an 
algebraic understanding of why the finite-T behaviour of bulk thermodynamic properties 
is mirrored in their /j,- dependence, and why meson masses decrease with \i even though /„. 
and — (qq) increase. The possibility of diquark condensation is canvassed. Its realisation is 
uncertain because it is contingent upon an assumption about the quark-quark scattering 
kernel that is demonstrably false in some applications; e.g., it predicts the existence of 
coloured diquarks in the strong interaction spectrum, which are not observed. 
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1. DYSON-SCHWINGER EQUATIONS 

The Dyson-Schwinger equations (DSEs) provide a Poincare invariant, continuum ap- 
proach to solving quantum field theories. There are many familiar examples, among them: 
the gap equation in superconductivity; and the Bethe-Salpeter equation (BSE) and co- 
variant Fadde'ev equation, which describe relativistic 2- and 3-body bound states. The 
DSEs are a system of coupled integral equations and a truncation is necessary to obtain a 
tractable problem. The simplest truncation scheme is a weak-coupling expansion, which 
generates every diagram in perturbation theory. Hence, in the intelligent application of 
DSEs to QCD, there is always a tight constraint on the ultraviolet behaviour. That is 
crucial in extrapolating into the infrared, and in developing uniformly valid, efficacious, 
symmetry-preserving truncations. 

The task of development is not a purely numerical one, and neither is it always obviously 
systematic. For some, this last point diminishes the appeal of the approach. However, 
with growing community involvement and interest, the qualitatively robust results and 
intuitive understanding that the DSEs can provide is becoming clear. Indeed, those 
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Figure 1. Q(k 2 )/k 2 from a solution [ ||] of the gluon DSE (dash-dot line) compared with 
the one-loop perturbative result (dashed line) and a fit (solid line) obtained following the 
method of Ref. [ £§]; i.e., requiring that the gluon propagator lead, via the quark DSE, to 
a good description of a range of hadron observables. 



familiar with the application of DSEs in the late- 70s and early-80s might be surprised 
with the progress that has been made. 

The DSEs have been used extensively [ [IJ in developing an understanding of confinement 
and dynamical chiral symmetry breaking (DCSB), and their wide application [ to the 
description of hadron properties in terms of their quark and gluon constituents is built 
upon that success. In understanding and unifying these phenomena, the DSEs point to 
the key role played by the necessary, momentum- dependent dressing of the elementary 
propagators and vertices in QCD. 

1.1. Gluon propagator 

Important in the application of DSEs is the gluon propagator, which at T = has the 
form: 

g ^ Ak) = ( 4 ^)W ff( * ):= _*L_, (1) 

where U(k 2 ) is the vacuum polarisation that contains all the dynamical information about 
gluon propagation. Studies of the gluon DSE have been reported by many authors [0J with 
the conclusion that, if the ghost-loop is not significant, then the charge-antiscreening 3- 
gluon vertex dominates and, relative to the free gauge boson propagator, the dressed gluon 
propagator is significantly enhanced in the vicinity of k 2 = 0. The enhancement persists 
to k 2 ~ 1-2 GeV 2 , where a perturbative analysis becomes quantitatively reliable. In the 
neighbourhood of k 2 = the enhancement can be represented [ |J cis ci rc gularisation 
of l/k 4 as a distribution, and this behaviour is illustrated in Fig. |TJ. A dressed-gluon 
propagator of this type generates DCSB and confinement without fine-tuning. (We identify 
confinement as the absence of a Lehmann representation for coloured propagators, with 
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Figure 2. M(p 2 ) := B(p 2 )/A(p 2 ) obtained in solving the quark DSE. The solution of 
M 2 (p 2 ) = p 2 defines M E , the Euclidean constituent-quark mass. 

the obvious analogue for other coloured n-point functions.) 

1.2. Quark propagator 

At T = = in a covariant gauge the dressed-quark propagator has the form 



S(p) :-- 



1 



ry-p + Efjo) 27 • pA(p 2 ) + B(p 2 ) 
is the renormalised dressed-quark self energy, which satisfies 

E(p) = (Z 2 - 1) i 7 -p + W + / A /^(P - q)^S(q)K(q,p), 



(2) 



(3) 



where is the dressed-quark-gluon vertex, mf- is the current-quark mass, ( is the 

renormalisation point, and := J A c? 4 g/(27r) 4 represents mnemonically a translationally- 
invariant regularisation of the integral, with A the regularisation mass-scale. Using Q{k 2 ) 
similar to that depicted in Fig. [I] and current-quark masses corresponding to 
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(4) 



one obtains [ ^ the dressed-quark mass function depicted in Fig. ||. 

For light quarks (u, d and s) there are two distinct domains: perturbative and non- 
perturbative. In the perturbative domain the magnitude of M(p 2 ) is governed by the the 
current-quark mass. For p 2 < 1 GeV 2 the mass-function rises sharply. This is the nonper- 
turbative domain where the magnitude of M(p 2 ) is determined by the DCSB mechanism; 
i.e., the enhancement in the dressed-gluon propagator. This emphasises that DCSB is 
more than just a nonzero value of the quark condensate in the chiral limit! 

For a given flavour, the ratio Mf /m^ is a single, quantitative measure of the im- 

portance of the DCSB mechanism in modifying that quark's propagation characteristics. 
As illustrated in Eq. (|^), 
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Figure 3. First two orders in a systematic expansion [ [jj of the quark-antiquark scattering 
kernel. In this expansion, the propagators are dressed but the vertices are bare. 
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this ratio provides for a natural classification of quarks as either light or heavy. For 
light-quarks Cf is characteristically 10-100 while for heavy-quarks it is only 1-2. The 
values of Cf signal the existence of a characteristic DCSB mass-scale: M x . At p 2 > 
the propagation characteristics of a flavour with mi < M x are altered significantly by the 
DCSB mechanism, while for flavours with mi ^> M x it is irrelevant, and explicit chiral 
symmetry breaking dominates. It is apparent from Eq. (|5|) that M x ~ 0.2GeV~ Aqcd- 
This forms a basis for many simplifications in the study of heavy- meson observables. [ ||] 

1.3. Hadrons are bound states 

The properties of hadrons can be understood by studying covariant bound state equa- 
tions: the Bethe-Salpeter equation (BSE) for mesons and the covariant Fadde'ev equation 
for baryons. The mesons have been studied most extensively and their internal structure 
is described by a Bethe-Salpeter amplitude obtained as a solution of 

[V H {k-P% u = [ A [x H (q;P)]srK r t :(q,k;P), (6) 

Jq 

where XH(q;P) ■= S(q + )T H (q; P)S(q-); S(q) := diag(S u (q), S d (q), S s (q), . . .); q+ = q + 
rjp P, g_ = g — (1 — rjp) P, with P the total momentum of the bound state and observables 
independent of rjp; and r,. . . ,u represent colour-, Dirac- and flavour-matrix indices. 

In Eq. (|6|), K is the renormalised, fully-amputated, quark-antiquark scattering kernel 
and important in the successful application of DSEs is that it has a systematic skele- 
ton expansion in terms of the elementary, dressed-particle Schwinger functions; e.g., the 
dressed-quark and -gluon propagators. The particular expansion depicted in Fig.||, with 
its analogue for the kernel in the quark DSE, provides a means of constructing a kernel 
that, order-by-order in the number of vertices, ensures the preservation of vector and 
axial-vector Ward-Takahashi identities and hence Goldstone's theorem. 

1.4. Phenomenological applications 

Some of the many applications of DSEs to the calculation of hadron observables at 
T — = n are summarised in Refs. [ |2| so here we make only two brief observations. Using 
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Figure 4. ?/;-meson electroproduction cross section: solid line; the dashed line is the 



p- meson cross section for comparison. The data are from Refs. [ 10 



the DSEs one obtains a mass formula for pseudoscalar mesons that unifies the small and 
large current-quark mass regimes. For mesons composed of quarks with small current- 
quark masses the formula yields [ || the Gell-Mann-Oakes-Renner relation. However, 
when the current-quark mass of one or both of the constituents is large, the formula 
predicts [ |§ that the meson mass increases linearly with the current-quark mass. In 
model studies [ |5| the linear increase is dominant for masses as low as that of the s-quark. 
The approach has also been extensively employed in the study of scattering processes, and 
one example is the calculation [ [| of the cross section for the diffractive electroproduction 
of vector mesons. This application makes the striking prediction, confirmed by recent 
data, that although two-orders of magnitude smaller than the p-meson cross section at 
the photoproduction point, the if)- meson cross section is equal to that of the p-meson at 
Q 2 = 15 GeV 2 , Fig.|. 



2. FINITE TEMPERATURE AND CHEMICAL POTENTIAL 

The contemporary application of DSEs at finite temperature and chemical potential is 
a straightforward extension of the T = = \i studies. The direct approach is to develop 
a finite-T extension of Ansatze for the dressed-gluon propagator and then solve the quark 
DSE. Having the dressed-quark and -gluon propagators, the response of bound states to 
increases in T and \x can be calculated. As a nonperturbative approach that allows the 
simultaneous study of confinement and DCSB, the DSEs have a significant overlap with 
lattice simulations: each quantity that can be estimated using lattice simulations can also 
be calculated using DSEs. That means they can be used to check the lattice simulations, 
and importantly, that lattice simulations can be used to constrain their model-dependent 
aspects. Once agreement is obtained on the common domain, the DSEs can be used to 
explore phenomena presently inaccessible to lattice simulations. 
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The renormalised dressed-quark propagator at finite- (T, /x) has the form 
S^(p, u k ) = ij -pA(p, u k ) + Z74 tUkC(p, uj k ) + B(p, u k ) (7) 
where tu k := u k + ip, with u k = (2k + 1) ttT, /c G Z, and satisfies the DSE 
S -1 (p, u k ) = Z£ ■ p + Z 2 (274 + m bm ) + £'(p, o; fe ) ; (8) 
where the regularised self energy is 

£'(p, tDfc) = if ■ pT,' A (p, u k ) + Z74 cDfe £' c (p, cD fc ) + £' B (p, u k ) , (9) 
S^(p,cj fc ) = jf ^g 2 D^(p- q,u} k - uJi)^tr [V^^S(q,uJi)T^(q,ui;p,uj k )} , (10) 

J& := T E^-oo / A (0 and 7^ := -(Zf/tfiry • p, P B := Z 1; P c := -(Zx/w^ry*. 
The complex scalar functions: A(p,a)fc), B(p,u k ) and C(p,uj k ) satisfy: .F(p, oj k )* = 
J r (p,U- k -i) , T = A, B,C, and although not explicitly indicated they are functions only 
of \p\ 2 and Cj k . The dependence of these functions on their arguments has important con- 
sequences in QCD. It can provide a mechanism for quark confinement and is the reason 
why bulk thermodynamic quantities, such as the pressure and entropy, approach their 
ultrarelativistic limits slowly 

2.1. Chiral phase transition at p — 

One order parameter for the chiral symmetry restoration transition is the quark con- 
densate, defined via the renormalised dressed-quark propagator: 



(qq)° c = Jim Z,(( 2 , A 2 ) N c T tr Dirac [S^q)] . (11) 

A-+oo Jl,q 



Since tr Dirac [Sm=o(o)] oc B(p,u k ), a simpler and equivalent order parameter is 
X(t,h):=ReB (p = 0,Cu ); t := - - 1 , h := — , (12) 

which makes it clear that the zeroth Matsubara mode determines the character of the 
chiral phase transition. (An order parameter for confinement, valid for both light- and 
heavy-quarks, was introduced in Ref. [ DJ. It is a single, quantitative measure of whether 
or not a Schwinger function has a Lehmann representation, and has been used [ 12] to 
striking effect in QED 3 .) 

The chiral transition is completely characterised by two critical exponents: (3, 5, which 
can be extracted from the behaviour of the chiral and thermal susceptibilities: 



, IN dX(t,h) 



dh 



. dX(t,h) 
Xt(t,h) : = 



Ot 



(13) 



h 



As a function of T the susceptibilities have a peak, which defines the pseudocritical points: 
tp C , t*, and because the correlation length is infinite at a second order transition it follows 
that 

Xl° := Xh(t h pc , h) oc hT* , z h := 1 - i , X T ■= Xt^ h) oc h~ z < , * := -1 (1 - 0) . (14) 
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Figure 5. Behaviour of the susceptibilities in the DSE model of two- light-flavour QCD 
illustrated in Fig. p],[ [14]] which is a minor modification of that in Ref. [ |fl: u; — >■ 0.6 GeV, 
that provides a slightly better description of tt and K properties. The critical temperature 
is Tf = 152 MeV and the critical exponents take mean-field values: Zh = 0.67, z t = 0.33; 
i.e., (3 = 0.5, 5 = 3.0, as might be anticipated because at long-range the interaction 
~ const, in configuration space. The other curves illustrate the slopes characterising 
mean-field (labelled IR) and 0(4) critical exponents. 



Hence the critical exponents can be obtained by plotting log 10 x pc vs - l°gio K as depicted 
in Fig.|5|. This figure makes one thing abundantly clear: very small values of the current- 
quark mass must be used to obtain accurate values of the critical exponents. If one 
retains only those values 5 x 10 -4 < h < 5 x 10~ 2 then an apparently good linear fit 
yields: Zh = 0.78, z t = 0.40, or (3 = 0.36, 6 = 4.5, which are quite close to the values 
of the 0(4) model. The small values of the current-quark mass we require to accurately 
estimate the critical exponents are inaccessible in contemporary lattice simulations. 

2.2. T and fi nonzero 

This is a difficult problem and the most complete study [ [TJ|] to date employs a simple 
Ansatz for the dressed-gluon propagator that exhibits only the infrared enhancement 
suggested by Ref. [ ||] : 

g 2 D,M n fc ) = fo* - r^yf^) ^ y 4o ^ ' (15) 

with Q k = 2kirT and rj a mass-scale. It is an infrared-dominant model and does not 
represent well the behaviour of D^ip^VLk) away from |p| 2 + VL 2 k « 0. Consequently some 
model-dependent artefacts arise. However, there is significant merit in its simplicity and, 
since the artefacts are easily identified, the model remains useful as a means of elucidating 
many of the qualitative features of more sophisticated Ansatze. 

Using Eq. fll5|) and the rainbow truncation [T^ — > 7 M in Eq. fllODl, the quark DSE is [ |7j 

S~\p,uj k ) = S^faQk) + jTi 2 -y u S(p,uj k )>y u , (16) 
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and we see that the simplicity of the Ansatz allows the reduction of an integral equation 
to an algebraic equation. Its solution exhibits many of the qualitative features of more 
sophisticated models. 

In the chiral limit Eq. ([H]) has two qualitatively distinct solutions. The Nambu- 
Goldstone solution, with 




, otherwise ^ - I1 + J1 + -4^- I , otherwise 

where pi := (p,u)i), describes a phase of the model in which: 1) chiral symmetry is 
dynamically broken, because one has a nonzero quark mass-function, B(pk), in the absence 
of a current-quark mass; and 2) the dressed-quarks are confined, because the propagator 
described by these functions does not have a Lehmann representation. The alternative 
Wigner solution, for which 

B(p k )=0 , C(p k ) = 

describes a phase of the model with neither DCSB nor confinement. 

The relative stability of the different phases is measured by a (T, /independent vacuum 
pressure difference: B(T,fi). B(T,p) = defines the phase boundary, and the deconfine- 
ment and chiral symmetry restoration transitions are coincident. For fi = the transition 
is second order and the critical temperature is T c ° = 0. 159 77, which using the value of 
rj = 1.06 GeV obtained by fitting the tt and p masses corresponds to T c ° = 0.170 GeV. 
This is only 12% larger than the value reported in Sec. |2.1| , and the order of the transition 
is the same. For any /i 7^ the transition is first-order. For T = the critical chemical 
potential is /i° = 0.3 GeV, which is ~ 30% smaller than the result in Ref. [ 16 1. The 
discontinuity in the order parameter vanishes at /a = 0. 

The quark pressure is easily calculated and is depicted in Fig. ^ Confinement means 
that P q = in the confined domain. In the deconfined domain it approaches the ultrarel- 
ativistic, free particle limit, -PuR; at large values of T and \x but the approach is slow. For 
example, at T ~ 2T C °, or fi ~ 3//°, P q is only 0.5 Pur- This feature results from the persis- 
tence of momentum dependent modifications of the quark propagator into the deconfined 
domain, as evident with C ^ 1 in Eq. fllBp . The figure also highlights the "mirroring" of 
finite-T behaviour in the /^-dependence of the bulk thermodynamic quantities. 

The vacuum quark condensate is given by the simple expression [ [TJ 



OAT 'max M; / fZ \ 

- (qq) = V 3 ^ I dyy 2 Rel^\-y* + p-uf-2itiu, l j : (19) 

T = T/rj, p, = fx/rj; l max is the largest value of / for which cof max < (1/4) + /i 2 and this also 
specifies u lmax , A 2 = uf max - Q?, pi = (y,Q l + ip). At T = = //, (-(qq)) = r/ 3 /(80vr 2 ) = 
(0.11 1]) 3 . Obvious from Eq. ([19]) is that (—(qq)) decreases continuously to zero with T 
but increases with /x, up to a critical value of fi c (T) when it drops discontinuously to zero. 
This behaviour is driven by the combination /i 2 — uf and is exhibited by the calculated 
result depicted in Fig.0. In the chiral limit one also has the simple expression 
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Figure 6. The quark pressure, P q (T,p), normalised to the free, massless (or Ultra- 
Relativistic) result. 
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It too involves the combination p? — uif and without calculation Eq. ( p0|) indicates that 
f n will decrease with T and increase with /i, as exhibited by calculated result depicted 
in Fig. [7]. These results confirm those obtained [ [11], [16| with more sophisticated Ansatze, 
and in doing so provide for their algebraic elucidation: the behaviour is a consequence of 
the necessary momentum dependence of the dressed-quark self energy. 
The (T, /x)-response of meson masses is determined by the ladder BSE 



r M (Pfc; Pi) = — 7 Re i 7m s(p 



r M {Pi]Pe 



In 



(21) 



where Pe := (P,£li), with the bound state mass obtained by considering Pe=o- In this 
truncation the u>- and p-mesons are degenerate. The calculated mass of the ir- and p- 
mesons is depicted in Fig. ^. The behaviour is easily understood by again considering the 
chiral limit where the mass of the longitudinal component of the p-meson is 



Mi 



1 2 



4(p 2 



7T 2 T 2 ) 



(22) 



The characteristic combination \i — ir T appears again. It is responsible for the anti- 
correlation between the response of M p to T and its response to fi, and that M 2 rises 
linearly with T 2 for /i = 0, just like a gauge-boson Debye mass. The mass of the transverse 
component of the vector meson is insensitive to T and p. 

In a 2-flavour, free-quark gas at T = the baryon number density is ps = 2/i 3 / (37r 2 ) , by 
which gauge nuclear matter density, po = 0.16 fm~ 3 , corresponds to \i = fio := 260MeV = 
0.245 rj. At this chemical potential the algebraic model yields 



M p {p ) « 0.75M p (/i = 0) , M^Guo) « O.85M (p = 0) 



(23) 
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Figure 7. — (qq) , f n and function of ji for a range of values of T. is only weakly 

sensitive to changes in /i and T because of mutual cancellations between the effects of 
these intensive variables on —(qq) and f n . That is a result of Goldstone's theorem; i.e., 
the preservation of the axial- vector Ward-Takahashi identity. 



where Afy(// = 0) = 1.02 GeV for m s = 180 MeV. The study of Ref. [ indicates that a 



more realistic representation of the ultraviolet behaviour of the dressed-gluon propagator 
expands the horizontal scale in Fig. [8], with the critical chemical potential increased by 
25%. This suggests that a better estimate is obtained by evaluating the mass at fi' = 
0.20 rj, which yields 



M p (fi' ) « 0.85M p (/x = 0) , MM) « 0.90M^(/i = 0) ; 



(24) 



a small, quantitative modification. The difference between Eqs. ( [Z3"D and (pi) is a measure 
of the theoretical uncertainty in the estimates in each case. Pursuing this suggestion 
further, fi 2 = /Uo\/2, corresponds to 2p , at which point M w = M p w 0.72 M p (/i = 0) and 
w 0.85 M^(fi = 0), while at the T = critical chemical potential, which corresponds 

m 



to approximately 3p in Ref. 



M,„ 



0.65 M p (ji = 0) and w 0.80 M^(p = 0). 



These are the maximum possible reductions in the meson masses. 

This simple model preserves the momentum-dependence of gluon and quark dressing, 
which is an important qualitative feature of more sophisticated studies. Its simplicity 
means that many of the consequences of that dressing can be understood algebraically. 
For example, it elucidates the origin of an anticorrelation, found for a range of quantities, 
between their response to increasing T and that to increasing \i. That makes clear why 
the transition to a QGP is second order with increasing T and first order with [i. Further 
it provides an algebraic explanation of why the (T, /^-dependence of ( — (qq)) and f n must 
be opposite to that of M p . 

3. DIQUARKS 

It is not implausible that diquark (quark-quark) correlations play an important role in 
the strong interaction. To see why, consider the class of models that can be characterised 
by an effective interaction of the form: 
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Figure 8. M p = and m n as a function of ft for T — 0, 0.1. On the scale of this figure, 
is insensitive to this variation of T. The current-quark mass is m = 0.011 rj, which for 
r] = 1.06 GeV yields M p = 770 MeV and = 140 MeV at T = = fx. 

/\ a \ a 
d 4 x d 4 y q{x)-f il —q(x) g 2 D flu (x - y) q{y)lv—q{y) ■ (25) 

Field theories defined by such an interaction admit [ [18] a meson-diquark bosonisation 
that at tree-level predicts a mass for both mesons and diquarks. The procedure corre- 
sponds closely to solving the Bethe-Salpeter equations obtained by combining the ladder 
truncation of the quark-antiquark and quark-quark scattering kernels, depicted in Fig.|3|, 
with the rainbow truncation of the quark DSE. Any form of g 2 D fiU (x — y) that is accu- 
rately able to reproduce the mass of the ir and p mesons in this truncation will predict [ 
[HJ the existence of stable, colour-antitriplet diquark bound states with masses: 

mli Rj 740 MeV, m# « 950 MeV, m\ A w mf- w 1500 MeV. (26) 

This points to a possibly important attraction in the quark-quark channel. 

However, it also presents a problem, of course, because asymptotic diquark states are 
not observed. The defect is in the truncation. [ If one proceeds beyond the ladder- 
truncation in the quark-quark scattering kernel; e.g., including the other diagrams de- 
picted in Fig. |3], one uncovers a repulsive contribution from the crossed-box diagram that 
eliminates the pole in the quark-quark scattering matrix; i.e., the full kernel does not 
support spurious diquark bound states. (The same procedure applied to SU(iV c = 2)- 
QCD predicts [ ^(J that the mesons and diquarks are degenerate in that theory.) It is 
clear therefore that analyses of diquark correlations based solely on an effective interac- 
tion of the type in Eq. fl2"5| ) can yield erroneous results, and the suggestion of diquark 
condensation may be amongst them. 

The consequences of a putative diquark condensate are easy to elucidate. It is adequate 
to consider a simple separable model [ |2"]| 

g 2 D pu (p- q, LJ k - u{) := V D g(p 2 + u 2 k ) g{q 2 + uf) , g{k 2 ) := exp(-k 2 /A 2 ) , (27) 
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Figure 9. Quark and diquark mass gaps: A m , A qq . A qq = at T = 0.15A, A = 0.71 GeV. 



which preserves Goldstone's theorem and ensures the absence of a Lehmann representation 
for the quark propagator. This Ansatz has two parameters that can be fixed by requiring 
a good description of pion properties at T = = /i; e.g., D = 106/A 2 , A = 0.71 GeV 
yield = 93MeV and — (qq) = (238 MeV) 3 . A feature of the model is that it admits 
a semi-algebraic analysis and at fi = 0, with these parameter values, it has second-order 
deconfinement and chiral symmetry restoration transitions at T c > 130 MeV, while for 
T = these transitions are first-order at ii c = 330 MeV neglecting the effect of a diquark 
condensate. 

A meson-diquark bosonisation yields a tree-level auxiliary-field effective action whose 
extremum is determined by a pair of coupled equations for the quark and diquark mass 
gaps, A m and A, 



"qq- 



A 
A 



d 3 p g 2 (pf) 



l=— oo ' 

oo 



qq 



l=— oo ' 



(2vr)3 e q {v x 
d 
(2 



£q(Pl) 



1 



0J 



f + (^)+A*) 2 + 4r(ft 



> -fx} 



,(28) 
(29) 



with 



4M) 



(30) 



The first equation arises from the confined-quark contribution to the effective action, 
while the second is the contribution from unconfined diquarks. The preliminary result of 
Ref. [ [21| for the (T, /^-dependence of A m and A qq is depicted in Fig. |9|. The indications 
are that at T = there is a first order transition from a A m ^ phase to a A qq ^ 
phase at p, ~ 280 MeV. This behaviour persists until T ~ 100 MeV when only A m ^ is 
possible. There appears to be a tricritical point in the (T, /i)-plane. 

Clearly the formation of a diquark condensate is a simple and direct consequence of as- 
suming that attraction dominates in the quark-quark scattering matrix. However, equally 
clearly, the validity of that assumption requires further consideration, with a first subject 
being the effect of the often-ignored repulsive terms in the quark-quark scattering kernel. 
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4. CLOSING REMARKS 

The DSEs are an efficacious tool for studying the strong interaction. In this application 
they expose the qualitative importance and quantitative effect of the nonperturbative 
dressing of propagators and vertices in QCD. The modelling involved is based on this 
observation, and it is necessary because of the need to make truncations. Questions will 
always be asked about the fidelity of that modelling, and the answers will come as more 
is learnt about about the constraints that Ward and Slavnov- Taylor identities in the 
theory can provide. That approach has already been particularly fruitful in QED,[ |22[ 
and in the development of a systematic truncation procedure for the kernel of the quark 
DSE and meson BSE.[ [7], [21J In the meantime the judicious application of DSEs will 
continue to provide a flexible, intuition-building framework for the prediction, correlation 
and validation of observables. 
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